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Abstract 



o 

We study the effect of a cut-off on the speed of pulled fronts of the one dimensional 
reaction diffusion equation. We prove rigorous upper and lower bounds on the speed in 
terms of the cut-off parameter e. From these bounds we estimate the range of validity 
of the Brunet-Derrida formula for a general class of reaction terms. 



1 Introduction 

The reaction diffusion equation 



u t = u xx + f(u) 



is one of the simplest models which shows how a small perturbation to an unstable state 
develops into a moving front joining a stable to an unstable state. The reaction term f(u) 
adopts different expressions depending on the physical problem under consideration. One of 
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the most studied cases, is the Fisher reaction term [13] f(u) = u(l — u) for which the asymp- 
totic speed of the propagating front is c = 2, a value determined from linear considerations. 
A more general case was studied by Kolmogorov, Petrovskii and Piscounov (KPP)[T!] who 
showed that for all reaction terms which satisfy the so called KPP condition 

f(u) > 0, /(0) = /(l) = 0, f(u) < f(0)u (2) 

the asymptotic speed of the front joining the stable u = 1 point to the unstable u = point 
is given by 

ckpp = 2 v / 7 7 (0)- 

The evolution of localized initial conditions to the front of minimal speed was established 
in [lj for general reaction terms. Recent work has dealt with effects not included in the 
classical reaction diffusion equation ([1]), namely the effects of noise and of the finiteness in 
the number N of diffusive particles. It was suggested by Brunet and Derrida [10] that such 
effects can be simulated by introducing a cut-off in the reaction term. In the case of noise the 
cut-off parameter measures the amplitude of the noise while in the case of finite number of 
N diffusing particles the cut-off parameter e — l/N. They presented numerical evidence to 
support their conjecture. By means of an asymptotic matching Brunet and Derrida showed 
that for a reaction term f(u) = u(l — u 2 ) a small cut-off changes the speed of the front to 

7T 2 

c~2-- -. (3) 

(log e) 2 

In recent work it has been show that the Brunet-Derrida formula for the speed is correct 
to (9((loge)~ 3 ) for a wider class of pulled reaction terms and cut-off functions. [7J [11]. A 
completely different behavior is found when a cut-off is applied to a bistable reaction term 
or to a pushed front, in these two cases the cut-off changes the speed by an amount which 
has a power law dependence of the cut-off parameter [6j [12] . The validity of representing the 
finiteness in the number of particles in the diffusion process by a reaction diffusion equation 
with a cutoff, and the effect of noise in the reaction diffusion equation with a cut-off was 
proved rigorously in [S] and [T5], respectively. 

The purpose of this work is to prove rigorous upper and lower bounds for the speed of 
fronts for reaction terms of the form f(u)Q(u — e) where / satisfies the KPP condition Eq. 
(J2J) and G is the step function. The results obtained are valid for all e, in the limit of e — > 
the upper and lower limits coincide and are the Brunet-Derrida value. 

The paper is organized as follows: in Section 2 we set up the problem and state our 
main result (Theorem 12. ip . In Section 3 we consider a relaxed variational principle which is 
crucial in determining upper bounds on the speed of propagation of fronts. In particular we 
prove existence of a unique minimizer for the relaxed variational problem. In Section 4, we 
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give an explicit expression for the minimizer of the relaxed problem. At the same time we 
characterize (in closed form) the value of the minimum of the relaxed problem, which allows 
us to give an explicit upper bound for the speed of propagation of fronts for (pQ). In Section 
5, we prove our main result (i.e., Theorem 12. ip and, in particular we provide error bounds 
for c. 



2 Statement of the problem 

The bounds will be obtained starting from the integral variational principle for the speed of 
the fronts. It was shown in previous work that the speed of the front of the reaction diffusion 
equation [TJ with arbitrary reaction term f(u) obeys the variational principle [21 IS] 

n F(l)/3 + f?F(u(a))/a*ds 
c =sup2 — — g 5 (4) 

u(s) J (du/ds) ds 

where sq = 1/ ' g(u = 1) is an arbitrary parameter, 

F(u)= [ U f(q)dq. 
Jo 

and the supremum is taken over positive increasing functions u(s) such that u(0) = 0, 
u(sq) = 1 and for which all the integrals in fl3J) are finite. 

We shall now be interested in reaction terms f(u) with a cut-off e of the form 



if < u < e 

u - N(u) if e < u < 1, 



where N(u), the nonlinearity, is such that N(0) = N'(0) = 0. Moreover we will assume that 
f(u) satisfies the KPP criterion Eq. ([2~]) , which implies N(u) > 0. 
Our main result is the following theorem. 

2.1 THEOREM. Consider the reaction diffusion equation (TJP where the reaction profile 
satisfies (TJ|). Moreover, assume N(u) < B[u — e) 1+v , for e < u < 1, where rj > 0. Then, the 
speed of propagation of fronts of the reaction diffusion equation (QP, c, satisfies, 

^ 2 - c ^ o (»)' (5 » 

Here, cl is given explicitly by 

c L = 2sin0*, (6) 



where 0* is the first positive solution of the equation, 

1 

2 



tan^* = i| loge|. (7) 
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In particular, for e — > 0, we have, 

_ 2 7T 2 / 1 \ 

L | log e| 2 \ | log e | 2 / 

Remarks: 

i) The speed cl is the precise speed of propagation of fronts for a profile defined piecewise 
as fh = 0, for < u < e, fiiu) — u, for e < u < 1, and fiiX) = (illustrated in the figure 
below) . 

1 

0,8 
0,6 
0,4 
0,2 


£ 1 




Figure 1: the profile fa corresponding to the relaxed variational problem 

ii) Although in principle the interest is focused on small values of the parameter e, our 
expression for = 2 sine/)* is valid for any < e < 1. In fact, one can also consider the 
interesting case e — > 1. In that case, the profile fi is peaked around u — 1, which is the 
typical situation that arises in the propagation of flames (first studied in p2]). For the case 
of profiles f(u) peaked around u — 1, the speed of fronts is approximately given by 




czfk = a/2 J f(u)du, 

(see, [T7]; it turns out that this expression czfk for the speed of the travelling fronts is 
actually a lower bound to the actual speed c, see [SI H]), Using the ZFK expression for the 
profile fi(u), one has, 

czfk = « v^l-e), 

as e — > 1. On the other hand, as e — > 1, | loge| = | log(l + (e — 1)| ~ 1 — e, approaches zero. 
Using ([7]), we see that also tan0* w in this case, and we have shi(/>* w tan0* w 0*. Hence, 
from ([7]), 
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and thus, 

c = 2 sin(0*) « 20* = ^2(1 - 
which coincides with the ZFK value. 



3 Relaxed Problem 

In this section we prove the existence of a maximizer for the functional 

i=d + P° ^p±ds 

so JO s 



This functional corresponds to (j3J) where the reaction term has been replaced by f^. It has 
to be maximized over all positive increasing functions u(s) suject to the conditions u(0) = 
and u(s ) = 1. 

Consider the functional 

for functions u : E + — >■ [0, 1] with tt(0) = 0, u increasing and lim^oo u(s) = 1, i.e., we do 
not reqire that u(s) assumes the value 1 at some point s . We denote this set of functions 
by C. For any positive, increasing function u with u(0) = and u(sq) = 1 for some sq we 
set v(s) := u(s) for s < s and v(s) = 1 for s > s . The function t>(s) is in C and a simple 
computation shows that 

F{u) = 2G(v) . (10) 

Thus, the supremum of J-{u) over all functions u and all values of s is not larger than the 
supremum of 2Q(u). 

3.1 LEMMA (ZFK bound). The functional Q is bounded above, in fact 



3.2 REMARK. Note that as e gets close to 1 the right side of ( 07]) fends to ^ anc ^ ^ence 
i/ie supremum of yjF{u) which corresponds to the wave speed, is less than y2{\ — e) as 
e -> 1. 

Proof. Fix any function w in C with both, the numerator and the denominator of Q finite. 
Since the derivative of u is square integrable, the function u is Holder continuous. In fact by 
an elementary estimate 

\u{s) - u(s')\ < y/\^7\y/\\ U '\\ 2 . (12) 
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Using the scaling invariance of Q we may assume that u(s) = e. The numerator M(u) of 
Q{u) can be written as 

f°° M 2 _ 2 roo 2 roc i 

N{u) = / ^^ds = ^ds-e = 2 —ds . (13) 

Je s Je s Je s 



Now by Schwarz's inequality 



M{u) < 2( ^ds) 1/2 ( / n' 2 ^) 1 / 2 (14) 



s 2 



or 

1 M(uf f°° „ , 

7 777 \ ^ < / u' 2 ds . 15 

Since u(0) = 0, w(e) = e and since the numerator does not depend on u on this interval [0, e] 
we find that the denominator of Q{u) is smallest when u{s) = s and hence it is bounded 
below by 

POD 

s+ u' 2 ds . (16) 



Thus 

A/-( M ) 2 + £ A^) n 7 ^ 

^ H " Vcu) 2 + 4^( M ) + 4. 2 (17) 

Note that the right side as a function of N(u) is increasing. Further, since 

N(u) < ] —^- (18) 
the stated estimate is established. □ 
3.3 LEMMA. There exists a function u G C such that 

G(u) = sup G(v) =: M . (19) 

vec 

Moreover, if we normalize u so that u{e) = e then on the interval (0, e), u{s) = s. 

Proof. Let u n be a maximizing sequence, i.e., 

Q{u n ) ->• M . (20) 

By scaling we can assume that u n (e) = e. By ffl8|) the numerator Af(u n ) is bounded. Hence 
the denominator is also bounded. Thus 

oo 

u' r 2 ds <C (21) 

for some constant C independent of n. By (fT2"|) the functions u n are uniformly continuous. 
Since the functons u n are uniformly bounded, by Arzela-Ascoli we can pass to a subsequence, 
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again denoted by u n which converges uniformly on any finite interval of [0, oo) to some 
function u. This function is in C since the pointwise limit of monotone functions is monotone. 
Since u n is bounded it follows from the dominated convergence theorem that 

lim / ^ds= ^ds , (22) 

and 

/*oo rco 

liminf / u'^ds > / u' 2 ds (23) 
n ^°° Jo Jo 

by the weak lower semicontinuity of the L 2 -norm. Thus 

M = lim Q{u n ) < g{u) (24) 

n— >oo 

and hence Q{u) = M. Note that both, the denominator and the numerator of Q(u) are finite. 
Assume that 

lim u{s) = a < 1 . (25) 

Then u/a G C and 

- > 6M ' (26) 

and hence a = 1. Thus u <E C and w is a maximizer. □ 

It remains to analyze the maximizers u of the functional Q{u). To this end we relax the 
functional once more. Let v be a function with < v(s) < l,v (0) = 0, and monotone except 
on some open subinterval of [0, oo). We denote this domain by V. Consider the function 

u v (s)= [ max{t/(t),0}cft . (27) 
Jo 

Clearly, this function is monotone and u v > v pointwise. Clearly, Af{u v ) > Af(v) and 
IKIh < \\v'h- Thus 

<3{v) < G{u v ) , (28) 
and the maximizer found before, maximizes the functional in the larger class T>. 

3.4 LEMMA, let u be a maximizer of the functional Q{u). Then 

so = inf{s : u(s) = 1} (29) 

is finite andu'(so) = 0. Moreover, on the interval (0, e], u(s) = s, on the interval (e, so] the 
function u(s) is of the form 

u (s) = v / iAcos(i v / 2/ M - 1) logs + 5) (30) 
for suitable constants A and 5. Finally, on (s , 00), u(s) = 1. 
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Proof. Assume that so = oo. By scaling we can assume that u(e) = e. Moreover since the 
numerator of Q does not depend on the function on the interval [0, e] the denominator is 
smallest by choosing u(s) = s in that interval. In particular we have that u(s) > e on [e, oo). 
Further u(s) < 1 for all s. Pick any smooth function / with compact support in (e, oo). 
Then for t small enough u + tf eD and 

j 2 f-2M I u'f = (31) 



i.e., 

ii 

— + 2Mu" = (32) 
s 2 

in the weak sense on the intervall (e, oo). Any weak solution of this equation is a linear 
combination of s a± where a± = |(1 ± a/1 — 2/M). By Lemma I3TT1 M < 2 and we get two 
complex conjugate roots. Hence the real solutions are 

u(s) = v / iAcos(^ v /2 / M - 1) logs + 5) , (33) 

where A and 5 are constants. The form of u(s) contradicts the monotonicity of u(s) and 
hence Sq < oo. This establishes also the form of u(s) on (e, sq\. To see that u'(sq) = pick 
any nonnegative smooth function /, whose support is in a close vicinity of sq. Clearly u — tf 
is in T> for t > and small. Then 

G{u - tf) < Q{u) (34) 

and hence 

- 2M J u'f > (35) 

On the interval (e, sq) the function u satisfies the equation f l3T]) and hence the right side of 
f )35|) can be rewritten as 



/ 



30 u CI 

-f-2Mu'f')ds+ / -fds (36) 

■'so S 



which, using integration by parts, yields 

< -2Mu'(s )f(s ) + [ \f(s)ds . (37) 

J S S 

Since u is increasing we have 

< 2Mu'(s ) < -J— [ \f(s)ds . (38) 

J [SO) J s s 

Assume further that / is a non-negative, smooth, compactly supported function with /(so) = 
1 and set f n (s) = f((s — s )n + s ). A simple calculation shows that as n tends to infinity 

/ - 2 f(s)ds « \ [ X > f(s)ds . (39) 
Hence u'(so) = 0. □ 
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4 The maximizer 

In this section we determine explicitly the optimizer, whose existence was established in the 
previous section. 

4.1 THEOREM. The unique maximizer is given by 

, s if < s < e 

u{s) = { - - (40) 

Ay/s cos(0(s)) if e < s < sq 

with 

A = ^T^ s ° = - ( 41 ) 
cos(0*J e 

and 

0( s ) = I C ot(0*)log(^)-0* . (42) 
Here 0* is the first positive solution of the equation 



0*tan(0*) = ~| log(e)| . (43) 

Moreover, we have 

M =: sup Q{v) = Q{u) = 2 sin 2 0,. (44) 
Proof. In a first step we show that 

s = \ ■ (45) 

We know by Lemma 13.31 and Lemma 13.41 that there exists a maximizer with the following 
properties 

u{e) = e u(s ) = 1 u'(so) = . (46) 
Since u'(s) is continuous and u(s) = s for s < e we also have 

u\e) = 1 . (47) 

Moreover, on the interval [e, Sq] the function u(s) is positive and increasing and has the form 



u(s) = VsAcos^y-^:-llogs + <yJ (4£ 

where M = Q(u), the maximal value of the functional. Note that by Lemma [3.11 M < 2. 
Since u(e) = e and u'{e) = 1 we have, using (J4"g|) 



: = .4(c)s(iJA_n og£ + 5 | ( 19) 
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-^="Vil- lsi,i (iVil- llog£+ v ■ (50) 

Similarly, from the fact that u(sq) = 1 and u'(sq) = we get from ( |4 



1 = v^Acos(iJ^ - Hog s + 5) (51) 



1 = ^° A V^" lsin( ^vF" llogSo + 5) ■ (52) 

Next we prove ()43p . i.e., we calculate = M. A straightforward calculation yields for 
the numerator 



00 \u(s) 2 - s 2 } 







s 2 



ds (53) 



y (log ^ + -fL=(mn{^ - 1 log s + 25)-sin( ] J^- 1 logs + 25) | (54) 

Likewise, for the denominator 

2/ M / (s) 2 rfs = 2e + ^-log^ (55) 







2M 



+ ( (sin( \/ ¥ 1 log s ° + 25) ~ sin{ \l Vi~ Uoge + 25) I (56) 



+ ^ |cos(J^:-llog S o + 2(5)-cos(J^:-lloge + 25) 1 . (57) 



The equation 



then reduces to 



X l " is) " S \s = m(2 I™ u'(s) 2 ds 



s 2 







= 2e-^J^-l[sin(jA-llog So + 25)-sin(jA-llog £ + 25) ) (58) 



^ 1 cos (\l^- U °SSo + 25)-cos(J^-l\oge + 25) J (59) 



Using f l49ll52|) together with the double angle formulas for cosine and sine one easily sees 
that the above equation reduces to 



e-T L )(s- 1 77)=0. ( 6 °) 



s J V 2 — M 



and hence fl4"5l) is proved. 
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The next step is to calculate M. Note that ( 15T1) and ( 1521) now read 



1 = 7? 4cos( 4vi|- llog£ + <5) (61) 



1 = ^Ti- ls ^Tr n ^ + S) - (62) 



from which we deduce that 



ta„(i^_ 1 log £ -5) = - 7 i=. (63) 



Likewise from ( 149]) and ( 150]) we obtain 



tan(iyA_iiog e + < f) = -_^L = . (64) 



Using the addition formula for the tangent function yields 



£-1 



tan (l/^ - llo S £ ) = -tan(y - l|loge|) = - - . (65) 



If we set 



and note that 



Wil~ 1|l0g£| ^ 



we learn that 



1 I log el 

tan 0» = —== = (68) 
.12 i 20* 

V M 1 

which yields (143]) . Since 0* > and M is the maximum of our functional, we have to choose 
0* to be the first positive solution of ( 115]) . In particular we have that 0* < n/2. 
It remains to determine 5 and A. Subtracting ( 163]) from ( 164]) we find that 

tan(2<5) = (69) 

and hence 5 = Nn/2 where iVeZ. Note that as s ranges from e to 1/e, the function 



1 / 2 

\ / TT — 1 log S 

2 V M 6 

varies from —0* to 0*. The function u(s) is positive and increasing and hence, if we choose 
the constant A positive, we find that 5 = 2nN where JVeZ. Hence we may choose 5 = 0. 
The function \\ — Hogs can be conveniently be written as 



2V ~k ~ llogS = \ coi( t ) * lo ^~ ~ ( t ) * ( 70 ) 
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and the condition that u(e) = e yields the value for the constant A stated in Theorem 
14. 1[ Finally, equation (jUJ) for the value of M follows immediately from the first equality in 

flEED. □ 



5 Error Estimates: Proof of Theorem 2.1 



In section 3 we have determined the exact value of the propagation speed, cl say, of initially 
localized disturbances obeying ([T]) for a linear profile with a cutoff. In fact, if the profile is 
given piecewise by f(u) = 0, for u < e, and f(u) — u for e < u < 1, we have shown that cl 
is given exactly by 

c L = 2sin0* (71) 
where 0* is the first positive solution of the equation 

0*tan0* = ^| loge|. (72) 

Solving (172"]) for 0* in power series on 1/| loge|, and replacing it in fITT]) one finds that 



c L = 2- 



7T 2 / 1 



loge| 2 \|loge 12 



where the leading two terms account precisely for the Brunet and Derrida result (i.e., Equa- 
tion |3] in the Introduction). 

Here, we would like to determine error bounds when the profile f(u) is a KPP profile 
with a cutoff, in other words, when the profile is given piecewise by f(u) = 0, for < u < e, 
and f(u) < u for e < u < 1. If we write f(u) = u — N(u), for e < u < 1, the KPP criterion 
amounts to requiring that N(u) > 0. For such a reaction profile, we have that 

F(u)= / f(q)dq, 
Jo 

is such that F{u) = for < u < e, whereas 

i r 
F{u) = _ (w 2 _ e2) _ y 

for e < w < 1. For a KPP profile iV(tt) > 0, thus, 

< G{u) (73) 

where 

G(n) = i(. 2 - e 2 ) + . 
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Hence, using ( 173|) in (J4]), and, taking into account (JHJ) and ( )44l) . we see that in general for 
a KPP profile with a cutoff, the speed of propagation of fronts for an initially localized 
disturbance of ([1]), say c, satisfies, 

c < c L . 

On the other hand, we can also use the variational principle embodied in (jlj) to obtain a 
lower bound on c. For that purpose we use as a trial function in (j4]) the minimizer u of the 
functional G. After some simple computations, we obtain, 

Cl ' ~ J? /e («W<fa 

Here, we will find estimates on the difference c 2 L — c 2 for profiles that satisfy the bound, 

< N{x) < B(x - e) 1+v (75) 
for e < x < 1, where 77 > 1. The denominator can be calculated in closed form as follows, 
D =j 1 J\u'{s)fds = e + j l J\u\s)fds 

= e + * . ^ sin ^ r : 2 sin(0» - tf dt 

4cos 2 </i*sm 2 (/)* cos</>, J—<f>* ! 

= e + , 2 / ■ ^ sin t ft (1 - cos(20, - 2t)) dt 

4 cos 2 0* sm 2 <A„ cos </>, J— 0* v v ™* '/ 

= e + ; 2/ ■ ^ sin l* (20* - (1/2) sin 40,) 

4 cos 2 0* sir 6* cos (A* V ^* \ I J r*J 



1 



' 4cos ,j (A* sin t 



(20* + sin 20,). (76) 



On the other hand, using ( !75|) in the numerator of (1741) . the properties of the trial function 
u(s) (in particular the fact that this function is increasing), we can estimate the numerator 

as, 

N< B ^ [ 1/€ (^ COS( ^ eV +, 'sin(0 -0) dS (77) 
~~ 2 cos 0, sin 0, J e \ cos 0, / * s 3 / 2 ' 

Using the fact that s = exp (20 tan 0,) (which follows from (|42p and (j4"5|) above) and that 

we can write i/e = exp ((loge)/2) = exp (— | loge|/2) = exp (— 0*tan0*), we have that 



y/se = exp [(0 — 0,) tan 1 
and also that 

ds 

2a/€ tan 0, exp ((0* — 0) tan 0, 



s 3/2 

Changing the variable of integration from s to in (177|) . making use of these last two 
expressions, we find, 



Be [+• ( cos0 \ 1+v 
N < — — — / exp [-(0* - 0) tan 0*] ■ e sin(0* - 0) exp (0, - 0) tan ( 

COS^ 0* Jm \ COS 0* / 



(7c 
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Finally making the change of variables — > a = 0* — we get, 



Be ^ W0*-o) V + V 
iv < — - — / exp — a~tan0J e sin cr exp ( a tan 0*) da. (79) 

COS 2 0* Jo \ COS 0* / 



Hence, from ( 1761) and ( J79l) . we have, 



iV cos 0* sin 0* . , 

— <45 y * . . J, 80 

D ~ 20, + sin(20*) v ; 

with, 

1=1 exp —a tan 0* e smcrexp (<Ttan0* j ) da. (81) 

Jo V cos 0* / 

When e — )• 0, we have from ( 1431) that 0* w 7r/2, sin0» « 1, sin(20*) w and cos0* = 
0(1/| loge|). Thus, in order to control the difference c| — c 2 , all we have to prove is that 

/<o(l/|loge|). 

We can estimate / from above by dropping the e inside the factor in the integral above. 
Moreover, we write cos(0* — a)/ cos0* = cos a + tan 0, sin a. Thus, we have, 



I < J = / (cos cr + tan 0*sincr) +n sin a exp (— ot? tan 0,). da. 
Jo 



(82) 



We now split the integral over a into two parts. We denote by J\ the integral between 
and a and by J2 the integral between a and 20, . The value of a will be conveniently chosen 
later on. We will first estimate J\. We use: i) exp(— arj tan (ft*) < 1 (since tan0* > 0); ii) 
cos a < 1 and iii) sin a < a < a to get 

f a a 2 
Ji< (1 + a t&n^) 1+v a da = — (1 + atan0*) 1+ \ (83) 
Jo 2 

and, using the convexity of x — > x 1+v (since 77 > 0), we have, 

Ji < Z^ 1 [a 2 + a 3+r? (tan0,) 1+ ^] . (84) 

On the other hand, in order to estimate J2, we use the fact that < cos cr, sin a < 1 in the 
interval [0,0*] (recall that 0* < vr/2). We also use that exp(— x) is decreasing, and we get 
at once 

■h < (1 + tan0*) 1+r? exp(— c*?7tan0*) 20*, (85) 
Pick any < r < 1, and then choose a to be, 

a = (tan0*)-( 2+ '" +r )/( 3+ ") (86) 
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The idea behind this choice is that it will make J\ = o ( 1 / 1 log e | ) , and at the same time it 
will make J2 of smaller order. Now, one can easily check that, for any < r < 1, 

2 1 > 1 + r. (87) 



3 + 7] 

Now, since for e small, tan0* > 1, it follows from (BID , our choice of a (i.e., equation (1861) ). 
and (|87D that, 

Ji < 2 r '(tan^)-( 1+r ). (88) 
Finally, using f H3|) in f l88|) . and the fact that r > 0, we get the desired estimate, 

J i = (tt-i) • ( 89 ) 



loge| 

Also, with our choice of a, (186|) . 

atan^ = (tan0*) (1 - r)/(3+,?) , 

where r > 1 and 77 > 0. Using fl85|) . we see that J2 is exponentially small as a function of 
1/1 loge| (the first factor grows polynomially as a function of tan0*, while the second factor 
is exponentially small). Summarizing, we have proven that 

1 



< c? - c 2 < o 



loge| 



6 Appendix 

For the sake of completeness, in this appendix we prove bounds on c 2 L — c 2 in terms of the 
parameter 77, for KPP profiles. These bounds allow us to show that as rj — > 00, c\ — c 2 — > 0. 
Consider, 

J— (cos cr + tan^* sin a) +v sinaexp (— arj tan*/)*) da. (90) 
Jo 

Denote by 

H == a tan 0* — log(cos cr + tan 0* sin a) , (91) 

and notice that 

rr - dH ■ (l + tan 2 0j 

if CT = — — = sin a- ■ — ; > 0. 

da (cos a + tan 0* sin a) 

Using (I9D|) and (19"T|) . we can write, 

J= (cos cr + tan 0* sin a) exp (—r/H) sin ada. (92) 



which can be rewriten as, 

-20 



1 r^* 2 

J = / (cos a + tan 0* sin a) exp (—r]H) H a da. (93) 

1 + tan (pi J 
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We have reamarked before that on the interval (0,20*), both < cos a, sin a < 1, thus 
(coscr + tan0* sincr) 2 < (1 + tan0*) 2 . Moreover, using that (1 + x) 2 / (1 + x 2 ) < 2 for x > 0, 
we can finally write, 



Jo 

Recalling that H a > 0, making the change of variables — >■ H, and computing H(0) = 
and H (20*) = 20* tan0* we get, 



Hence, J — > as rj — > oo. 
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